Topological Excitations of One-Dimensional Correlated Electron Systems 



M.I. Salkola 1 * 2 and J.R. Schrieffer 1 
1 NHMFL and Department of Physics, Florida State University, Tallahassee, Florida 32310 
2 Department of Physics, Stanford University, Stanford, California 94305 

(August 11, 1998) 

Properties of low-energy excitations in one-dimensional superconductors and density-wave systems 
are examined by the bosonization technique. In addition to the usual spin and charge quantum 
numbers, a new, independently measurable attribute is introduced to describe elementary, low- 
energy excitations. It can be defined as a number w which determines, in multiple of n, how many 
times the phase of the order parameter winds as an excitation is transposed from far left to far 
right. The winding number is zero for electrons and holes with conventional quantum numbers, but 
it acquires a nontrivial value w = 1 for neutral spin-i excitations and for spinless excitations with 
a unit electron charge. It may even be irrational, if the charge is irrational. Thus, these excitations 
are topological, and they can be viewed as composite particles made of spin or charge degrees of 
freedom and dressed by kinks in the order parameter. 



The concept of elementary excitations due to Landau 
plays a fundamental role in understanding many proper- 
ties of condensed-matter materials. It relies on the as- 
sumption that in the renormalization-group sense there 
exists a map onto an effective model with the same low- 
energy, long- wavelength physics but with few relevant de- 
grees of freedom left. The practical significance of this 
approach is that seemingly complex systems can be stud- 
ied and experimental predictions made even though true 
microscopic interactions may be strong. 

In spite of the success of the basic concept, it cannot 
be justified rigorously because of the break-down of the 
renormalizing procedure. Indeed, there are many inter- 
esting instances where the physical picture of weakly in- 
teracting elementary excitations with quantum numbers 
equal to those of bare ones becomes invalid. Examples of 
such cases are the quantum Hall effect and quasi-one- 
dimensional conductors Q where low-energy excitations 
have unusual spin-charge relations and are not continu- 
ously related to conventional electrons and holes. 

In this Note, we examine interacting one-dimensional 
conductors. This is a widely relevant and frequently 
studied problem [||J^] . We consider a particular situation 
in which an extra particle is added into a superconduc- 
tor. Yet, the consequences of this process illustrate that 
the elementary, low-energy excitations are more complex 
than one might have argued. Our most important result 
is that, in addition to the usual spin and charge quantum 
numbers, an excitation attaches itself a kink which can be 
quantified by another quantum number. It is defined as a 
number w which determines, in multiple of 7r, how many 
times the phase of the order parameter at a given point 
winds as the excitation is transposed from one to another 
end of the system. While the winding number w does not 
constitute an independent degree of freedom in the sense 
that it would be unrelated to the spin and the charge of 
the excitation, it does have physical attributes that make 



it measurable without any information about them. The 
winding number would be zero for an electron and a hole 
with conventional quantum numbers, if they were to ex- 
ist as elementary excitations, but it acquires a nontrivial 
value w = 1 for a neutral spin-^ excitation and for a spin- 
less excitation with a unit electron or hole charge. It may 
even be irrational, if the charge is irrational. Thus, these 
excitations are topological objects which can be viewed 
as composite particles made of spin or charge degrees of 
freedom and dressed by e i7 ™-phase kinks in the order pa- 
rameter. The winding number appears naturally in sys- 
tems which have continuous symmetry and whose ground 
states develop quasi-long-range order in one dimension 
and true long-range order in higher dimensions. In addi- 
tion to superconductors, kinks must then be introduced 
in spin- and charge-density-wave systems. Our conclu- 
sion complements the earlier observation || that soliton 
excitations have unusual fermion quantum numbers in 
charge-density-wave systems where massless fermions are 
coupled to a boson field. Here we demonstrate that also 
the converse is true: elementary charge and spin excita- 
tions always carry kinks in one-dimensional systems with 
quasi-long-range order. 

Specifically, consider the one-dimensional Hamiltonian 

H = -t 2j(4 +lo .CW + h.C.) + \ ^2 Vn-mPnPm, (1) 
na nm 

where c na is the fermion operator for an electron of 
spin a at site n, p n = c^cw is the electron num- 
ber density, V n is the electron-electron interaction, and 
t is the nearest- neighbor hopping matrix element. To 
address the low-energy and long- wavelength excitations, 
the energy spectrum is linearized at the Fermi energy and 
the fermionic degrees of freedom are expressed by slowly 
varying fields ip a ± (x) : 

<W>/S~ ^+(x n )e- ikFX " +i_(. n )e ,tFl ", (2) 
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where kp is the Fermi wavevector, a is the lattice spac- 
ing, and x n = na. Subsequently, the free part of the 
Hamiltonian, H = Hq + Hj, becomes 



Hq = v F Jdx [V4(x-)p?Acr(x-) 



^(x + )pVv(x+)], (3) 



where p — —ifid x is the momentum operator and 
vp = (2ta/H) sinkpa is the Fermi velocity (hereafter 
a sum over repeated indices is implied unless other- 
wise noted). The left (+) and right (— ) moving elec- 
trons with spin a are labeled according to their argu- 
ments, denoting both the space and time coordinate, 
x± = (x,t)±; namely, Vv(x±) = ipa±{x,t). In the case 
of free fermions, the Heisenberg equations of motion, 
ihd t ^pa(^±) = [Vv(x±), H], lead to the field operators 
which are functions of variables x ± vpt, where ± re- 
fer to the left- and right-moving electrons, respectively: 

-0<t(x ± ) = lp a (x ± Vpt). 

It is convenient to define left- and right-moving cur- 
rents as J aa '(x±) — :^pl(x±)ip cr r(x±):. The colons de- 
note normal ordering with respect to the filled Fermi sea 
of the noninteracting system. In terms of these currents, 
the Hamiltonian describing the electron-electron interac- 
tions may be rewritten as 

Hi = Jdx (-Vi J CTp (x + )J jLlo .(x_) + V 2 J crC r(x + ) J w (x_) 

+iy 3 [e*^t(x.)^(x + ): :^(x_)^(x+): + h.c] 
+ iF 4 [J™(x + )J W( (x + ) + J ff(T (x_) J MM (x_)]) . (4) 

As usual, V\ is the backward-scattering constant, V% 
is the forward-scattering constant, V3 is the Umklapp- 
scattering constant, and V4 is another forward scatter- 
ing constant. Because V3 corresponds to the scatter- 
ing processes which violate momentum conservation by 
Akp, it is important only if Akp is equal to the recip- 
rocal lattice constant (the second-order commensurabil- 
ity). Note that Fi >3 = V(2k F ) and V 2A = V(0), where 
V{k) = J2n aV n e * kna - For simplicity, we consider a 
limit where the backward and Umklapp terms are ei- 
ther zero or scale to zero, so that the model is exactly 
solvable by bosonization. For example, far away from 
the second-order commensurability, the Umklapp pro- 
cesses are effectively turned off and remain so under the 
renormalization-group flow. As long as the backward and 
Umklapp scattering terms remain irrelevant in the sense 
of scaling, they can be neglected. Finally, the coupling 
constants T4 are allowed to be independent, as implied 
by many important interactions that are not depicted by 
the original Hamiltonian, Eq. (El). 

In Abelian bosonization |3|||f^], the left- and right- 
moving fermion operators are expressed in terms of bo- 
son fields, Vv(x±) ~ if CT :exp[q=i$ (T (x±)]: (no sum over 
repeated indices is implied in this paragraph), where K„ 
is the Klein phase-operator which establishes the correct 
anticommutation relations for different fields B. For 



instance, with this definition, the current operator be- 
comes J CTCT (x ± ) = (l/2^)5 x $ CT (x±). For Vy = U3 = 0, 
the Hamiltonian is quadratic in the boson fields 

Ho = Y,hv F f ^ ([a c <Mx+)] 2 + [5A(x-)] 2 ) , (5a) 

Hi = T,[ ^(m^wnwx-)] (5b) 

+ iu 4 {[a c <Mx+)] 2 + [a a $ M ( x _)] 2 }). 

In the absence of interactions, $ CT (x±) = $o-(x ± VFt). 
H is readily diagonalized by a set of unitary transforma- 
tions. First, define $ s (x±) = [$-|-(x±) - $j.(x±)]/\/2 
and $ c (x±) = [<&-[• (x±) + c f ) |(x±)]/v / 2- Second, de- 
fine (p s (x±) = $ s (x±) and y c (x±) = $ c (x±) cosh</) + 
$ c (x=f:) sinh (j), with g — Vijihvp + V4) and tanh20 = g. 
The Hamiltonian is transformed into 

H=J2 hv »[ % (I^( x +)] 2 + I^(x-)] 2 ) , (6) 

where v s — vp and v c — Dp/ cosh 29 are the spin and 
charge velocities. In the Heisenberg picture, </? s (x±) = 
ip s (x ± v s t) and <p c (x±) = (p c (x ± v c t). A useful relation 
in computing correlation functions of interacting fields is 
the two-point correlation, ((p u (x.±)ip„(0±)) = — log(a =p 
ix u ±), where x„± — x ± v v t (y — s,c). The ultraviolet 
cutoff, the Fermi length, is formally associated with the 
lattice spacing a. 

For g < 0, the superconducting (SC) instability is the 
most dominant one at zero temperature. Defining the 
singlet pairing field as 

A(x) = [^(x+)^ T (x_) - Vt(x+)V>|(x-)]/\/2, (7) 

the correlation function 

4°)(x-y) = (0|At(x)A(y)|0) (8) 

probes the ordering fluctuations of interest in the ground 
state |0). Below, we will always consider cases where x 
and y are measured at equal times t: x = (x, i) and 
y = (y,t). Using the identity, 



one can show that 



-a/3(Vj(x) ¥ >(y)> 



, (9) 



At(x)A(y) = G sc (x - y) cos A.:, (10) 

where the prefactor is 



Gsc(r) 



," / ft 2 \ Msc/2 



(11) 



The operators, A v = [<^„(x) - tp„(y)]/y/2 (y = s,c), 
are expressed in terms of new fields, y s (x) = ip s (x + ) + 
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¥> s (x_) and <y5 c (x) = <y9 c (x + ) — ip c (x^). The exponents 

are defined as fisc = 1 + »7 2 > A* = 1 + V\A — 9 2 > 
and 77 = [(1 + <?)/(l — ff)] 1 ^ 4 - In the ground state, 
Gg^x) = Gsc( x ) — thus, the asymptotic forms of the 
correlation function are Gg^(x) cx x~^ sc and G^(q = 
0,o;) tx oj~ qsc , where asc = 2 — /Use- The supercon- 
ducting correlation function is singular (asc > 0) at low 
energies and small momenta, for g < 0. 

The effect of excitations on pairing fluctuations can 
be examined by computing the correlation function in 
the state |* T ) = ip\{0)\0), where ^{(0) = [ip\{0+) + 

V , |(0-)]/\/2 adds an electron into the system at the origin 
and time t = 0. First, let the electron be far away from 
the points x and y where the pairing fields are measured, 
\r\ <C |-R|; here, r = x — y and R — (x + y)/2. Then, writ- 
ing A v = (r / y2)d x ip v (R) + C(r 3 ), an operator-product 
expansion can be developed for Eq. ( |l0| ) . The correlation 
function, 



GW(x,y) = <* T |At(x)A(y)|* T >, 



T sc\ 

simply becomes for t = 



G^(x,y)=l-1 



1 v 3 M /2 



COS 



(12) 



ar/2 
a 2 + R 2 ' 

(13) 



For given r, the influence of the electron on the correla- 
tion function decays as R~ A . This behavior is universal 
as the scaling exponent is independent of the strength 
of interactions. The result is consistent with the natural 
observation that a single electron has no effect on the 
bulk properties of a superconductor. 

Second, let x and y be arbitrary. The correlation func- 
tion can be written in the form 

G«(x,y) =G SC (r)F(x,y), (14) 

where F(x,y) = [F+(x, y) + F_(x, y)]/2 and 



^±(x,y) = cos[6 a (x s ±) - 8 a (y s ±)] 



(15) 



xe ±iw 1 [e a (x c± )-e a (y c± )]±iw 2 [9 a (x cT )-e a (y cl :)]_ 

We have defined 9 a (x) — arctan(x/a), w% = 77 cosh tfi, 
and u>2 = — r)svah<p (w\ > W2 > 0, for g < 0). Note 
that # a (£) — ► f sgn(x), as \x\/a — * 00. Initially (rj = 0), 
the function F reduces to F(x, y) = cos 2 [0 a (a;) — 6 a (y)]. 
For 1 7- 1 -c \R\, we recover Eq. (|13|). In the scaling limit 
where r, i? — > 00 with z = r/R < 1 fixed, the asymptotic 



behavior of Gg 1 ^ is 



Gsc(x,y) cx 



i? 



(16) 



where A/a = [2z/(l — z 2 )] 2 . Again, the scaling expo- 
nent describing the suppression of superconducting cor- 
relations near the electron is universal. 
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FIG. 1. Schematic illustration of propagating kinks in the 
order parameter A(x) tx F(x, 00) after a right-moving elec- 
tron has broken up into spin and charge excitations in a su- 
perconductor away from half filling; Vi = and V2 < 0. The 
left-moving charge component is not shown. Note that, in the 
interacting system, the spin velocity is larger than the charge 
velocity, v s > v c . 

In general, the order parameter A is obtained from 
(At(x)A(y)) — > |A| 2 , as r — > 00. Likewise, the influ- 
ence of an excitation on the superconductivity can be 
studied by letting x to probe the order parameter in the 
neighborhood of the excitation and y — > 00, even though 
strictly speaking A vanishes in one dimension. For t = 0, 
the solution is 



F(x,oo) = 1 - 



(17) 



which shows that the electron destroys the order param- 
eter in its vicinity || while preserving F > 0. The 
order-parameter suppression dies out as x~ 2 . However, 
in the interacting system, the electron decays instan- 
taneously into elementary excitations involving cither 
spin or charge degrees of freedom. The kinky nature 
of these excitations is revealed, when the system is al- 
lowed to evolve in time. For illustrative purposes, con- 
sider the ground state with one right-moving electron 
initially added at the origin, |* T ) = ?/){(0_)|0), so that 
F = F_ . The spin and the charge excitations carry kinks 
in the order parameter, and they both separately con- 
tribute to the total phase shift an equal amount of n, 
because w± + W2 = 1. Initially, the kinks overlap but 
then split as the excitations move apart. If only the spin 
excitation is located between the points x and y at a 
later time t, the correlation function is negative: G^ 
changes sign if either x or y crosses the position of the 
spin. This clearly shows that the state has an or- 
dinary kink, or an antiphase domain wall, in the order 
parameter precisely at the position of the spin fT(| . The 
winding number of a spin-^ particle is 1. In contrast, 
the splitting of the charge into right and left moving 
charge components, ql and q^, leads to two irrational 
kinks whose winding numbers are w\ and W2 ■ The wind- 
ing number of the kink is uniquely related to its charge 
— in units where the electron charge is one, the relation 
is particularly simple: — q\2- The total winding 

number, as well as the total charge, is conserved. Be- 
cause the spin and charge velocities differ, the spin and 
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charge kinks also propagate with the different velocities. 
This suggests a new and potentially attractive way to de- 
tect spin-charge separation in a superconductor by using 
the Josephson effect. It would also serve as evidence of 
irrational charge. Figure 1 describes schematically the 
effect of spin and charge on the superconducting order 
parameter. 

The above result is generalized in a straightforward 
manner for N electrons initially injected into the system 
at positions Rj = (i?i, 0) (i — 1, . . . , N): 

N 

4c ( x > y) = G sc (r)Y[F( X ~R l ,y-R l ). (18) 

8=1 

If an infinite number of electrons is randomly distributed 
with a mean distance £, the correlation function decays 
exponentially: 

G^(x,y)oc^y SC e -M^ (19) 

with r(t) = min[|r , |, (v s — w c )|t|]. Thus, for any nonzero 
concentration of injected electrons, n ex = aji 7^ 0, su- 
perconductivity is destroyed. 

One may also ask whether other quantities than the or- 
der parameter show any signatures of kinks. Indirect evi- 
dence of kinks could be looked for in the density of states 
and the conductivity, for example. They are probed in 
specific heat and optical absorption measurements, (i) 
That there are new states at low energies ("midgap" 
states) associated with the kinks is evident from the den- 
sity of states N(lo). In the absence of excitations, at zero 
temperature, N(oj) ~ w°, where a — |(l/yT — 9 2 ~ !)■ 
Enhanced superconducting fluctuations will always lead 
to a pseudo-gap in the density of states, because a > 0, 
for g < 0. If there are excitations present in the system, 
the density of states must correspondingly be modified 
at low energies such that 

N(fij) ex (to 2 + 7l 2 )"/ 2 + (co 2 + 7 |)"/ 2 , (20) 

where 71 = (2v c /£) [cosh 2 (j> — sin(^ cosh 20)] and 72 = 
(2v c /£)[l - cos(§ sinh20)]. Thus, nonzero N(lo -> 0) 
implies that kinks give rise to new states at low ener- 
gies, (ii) The conductivity, however, does not exhibit 
any unusual behavior due to kinks, because it is sensi- 
tive to the concentration of charges and the interactions 
between them — injected electrons behave undistinguish- 
ably from the existing particles, and the charge dynamics 
remains the same. If kinks and excitations were pinned, 
different kind of behavior would be expected. 

In general, a nonzero backward scattering term which 
scatters electrons of opposite spin across the Fermi sur- 
face in opposite directions must be included. As a result, 
left- and right-moving spin degrees of freedom are cou- 
pled so that the spin will also split into left- and right- 
moving components. While it appears that the spin by 



itself cannot have other values than integers and half in- 
tegers, irrational values are allowed for the average spin 
and the concomitant kink. 

The concept of a winding number applies equally to 
charge-density-wave and spin-density-wave instabilities 
(at 2k f) that are described by operators 

p(x) = ^(x+)W(x_)e i2fc ^ + h.c. (21a) 

(charge-density-wave) and 

S„(x) = ^(x+)r CT V p S v(x_K 2fe ^ + h.c. (21b) 

(spin-density- wave) ; r a are the three Pauli matrices (a = 
1,2,3). The only difference is that the functional forms 
of the irrational winding numbers w depend in a unique 
fashion on the nature of quasi-long-range order. 

In conclusion, elementary excitations in interacting 
one-dimensional conductors always carry kinks in the or- 
der parameter. The winding number characterizing a 
kink cannot be arbitrary but is determined by the spin 
and the charge of the excitation and by quasi-long-range 
order of the ground state — in other words, by the in- 
teractions. This is not an entirely unexpected result, be- 
cause electrons are interpreted in bosonization as solitons 
S . The novelty of our formulation is that the kink struc- 
ture of the excitations becomes observable, if the ground 
state develops quasi-long-range order. This implies a 
completely new conception of probing unusual quantum 
numbers of elementary excitations, which manifest them- 
selves through structural and dynamical modulations in 
the order parameter. For example, in a superconduc- 
tor, the Josephson effect could provide a method to mea- 
sure winding numbers associated with injected elemen- 
tary excitations. Using Hartree-Fock theory, we have 
already pointed out that spin excitations in supercon- 
ductors tend to form antiphase domain walls (kinks) in 
the order parameter [j]~T]| . An analogous observation con- 
cerning impurity spinons has been made in the context 
of one-dimensional Kondo lattices Thus, these re- 

sults corroborate the conclusion that kinks are generic 
excitations of superconductors. 
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